Meter-baseline tests of sterile neutrinos at Daya Bay 
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We explore the sensitivity of an experiment at the Daya Bay site, with a point radioactive source 
and a few meter basehne, to neutrino oscillations involving one or more eV mass sterile neutrinos. 
We find that within a year, the entire 3+2 and 1+3+1 parameter space preferred by global fits can 
be excluded at the 3a level, and if an oscillation signal is found, the 3+1 and 3+2 scenarios can be 
distinguished from each other at the 3a level provided one of the sterile neutrinos is lighter than 
0.5 eV. 



Introduction. The standard three neutrino (3;^) pic- 
ture has been successful in explaining most oscillation 
data. However, data from the Liquid Scintillator Neu- 
trino Detector (LSND) experiment ^Ij when interpreted 
as arising from — 9^ oscillations, indicate a devia- 
tion from the simple picture. The Mini-Booster Neu- 
trino Experiment (MiniBooNE) |51 provides supporting 
evidence for the LSND result that oscillations involving 
an eV mass sterile neutrino may be at work. Additional 
support may be found in an upward revision in the esti- 
mate of the reactor flux yield p]. The fact that short 
baseline (SBL) reactor neutrino experiments do not de- 
tect the 3% larger flux could be explained as a conse- 
quence of oscillations to sterile states. 

Popular scenarios that are consistent with the the rel- 
evant data have either one sterile neutrino, with a 3+1 
mass spectrum (such that the nearly degenerate triplet 
of mass eigenstates is lighter than the remaining state), 
or 2 sterile neutrinos [11 15]. The 5 neutrino (Bi/) case has 
2 viable spectra: a 3+2 spectrum in which the triplet is 
lighter than both sterile neutrinos, and a 1+3+1 spec- 
trum in which one sterile neutrino is lighter than the 
triplet and one is heavier. In all cases, the sterile neutri- 
nos mix little with the active neutrinos. 

Recently, it was suggested that a ten kilocurie scale 
-^''^Ce-^'^^Pe /3-decay source could be placed inside a large 
liquid scintillator detector to study eV sterile neutrino os- 
cillations on baselines of a few meters with 1.8-3.3 MeV 
neutrinos [6.. Distinct virtues of this technique are (1) 
that with a point-like source, an oscillation signature can 
be demonstrated as a function of both energy and base- 
line, (2) the short baseline may be easily adjustable, and 
(3) existing detectors can be utilized. Clear technical 
challenges are the feasibility of constructing such an in- 
tense radioactive source and of engineering suitable ultra- 
pure shielding of the source inside the detector. For a 
decisive measurement, Ref. (7). considered the possibil- 
ity of an experiment at the Daya Bay site with a 500 kCi 
source. The configuration of the 4 detectors in the Far 
Hall at Daya Bay makes it possible to place the source 
outside the detectors thus circumventing one of the tech- 
nical issues. 

In this Letter we show that the parameter space pre- 



ferred by global fits in the 3+1, 3+2 and 1+3+1 scenar- 
ios will be stringently tested by the proposed multi-meter 
baseline Pe disappearance measurement at Daya Bay. For 
sterile neutrino masses below 0.5 eV, such a measurement 
can even distinguish between the 3+1 and 3+2 scenarios 
at the 3(7 level. This enhanced sensitivity arises because 
knowledge of the fraction of the 1/4 and 1^5 mass eigen- 
states breaks the degeneracy in the sterile mixings to Ve 
and t'p, both of which are required to explain the anoma- 
lous SBL data. 

Sterile neutrino oscillations. For vacuum oscilla- 
tions of MeV neutrinos from a radioactive source, the 
(CP phase-independent) i/g and survival probability 
at distance L is 

Pee = 1 - 4 ^ \Uej |' Siu^ A,, , (1) 

where Aij — Srn^^L / lAE^) with 6m1^ — ml — m|. i,j 
denote the mass eigenstates and take values from 1 to the 
total number of neutrinos. Uei are elements of the mixing 
matrix. For the 3+1 spectrum, i5to43 ~ 6m'\2 — Srrvl^ ~ 
1 eV^ > (5m§2 ^ (5m§i ~ 2.4 x lO'^eV^ > Sm^^ ~ 
7.5 X lO-^eV^ Then, P^+^ = 1 - sin^ 261^ sin^ A41, with 
the definition, sinOs = Uei. 

In the 5u case, Eq. ([T]) includes a superposition of three 
oscillation frequencies corresponding to Snili,Sm^-^ and 
(5m54. We neglect mixing between the 2 sterile states, 
so that Sml^ does not contribute. Although the sterile 
neutrinos can mix with all three active neutrinos. Pee 
depends only on the four parameters, Srrili,Smli, \Ue4:\ 
and \Ue5\ via 

P!: = l-i{l-\Ue4\'-\Ue5n 

X (|[/e4|'sin2 A41 ^ \u^^f sin^ A^^) . (2) 

Since P^^ is insensitive to the signs of A41 and A51, 
Ve disappearance data cannot distinguish between the 
3+2 and 1+3+1 spectra for identical mixing matrix el- 
ements. (If we allowed for mixing between the sterile 
states, the right hand side would have an additional term, 

-4|C/e4p|f/e5|'sin2A54.) 

In the left panel of Fig. [l| we show the i/e survival 
probability for several 3+1 and 3+2 sample points. For 
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FIG. 1. Left: The energy-averaged Ve survival probability 
as a function of distance for 3+1 and 3+2 sample points. 
Ue4 ~ 0.16 (giving a ~ 10% oscillation amplitude), and in 
the 3+2 scenario, Ues is also 0.16. Right: Event distributions 
for the chosen radioactive source-detector configuration. The 
solid and dashed curves show the cases of no active-sterile 
oscillations fTI, and of oscillations with Sm^ — 1 eV^ and a 
10% oscillation amplitude, respectively. 



the sake of illustration, we have used somewhat large 
values of Uei and J7e5- The significant variation in the 
survival probabilities over the first few meters for dif- 
ferent {Smli,Smli) choices reveals the strength of the 
method. For all curves in Fig. [T] Pee is convolved with 
the energy spectrum from the radioactive source. 

Experimental set-up and procedure. The 500 kCi 
radioactive source at Daya Bay can be placed so that 
the 4 cylindrical detectors collect data with baselines 
from 1 to 8 meters. Several possible source locations have 
been studied, each giving a different spatial coverage of 
PeeiL). We choose "Point B" in the jargon of Ref. [7], 
which is located halfway between two of the detectors, 
and samples 2 principal baselines. It provides superior 
sensitivity for Sni^ ~ 1 eV^ with an oscillation length of 
about 6.5 meters. The no oscillation signal event rate 
is about 38,000 per year. Event distributions as a func- 
tion of baseline are shown in the right panel of Fig. [T] 
Depending on the energy window used, the reactor neu- 
trino background is expected to He between 22,000-32,000 
events per year. However, this large background can be 
controlled because its shape will be known. 

We take the detectors to be identical and adopt the 
following for our analysis [J: 
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where iVf ? is a simulated dataset and Nf^ is the theoret- 
ical expectation for a given set of oscillation parameters, 
and i and j run over position and visible energy bins, 
respectively, (jg — 0.01 and ar — 0.01 are the normaliza- 
tion uncertainties in the signal and reactor background 
fluxes, respectively, and cTf, = 0.02 is the bin-to-bin un- 
certainty, as and are nuisance parameters that are 
allowed to float. N^^ and A^*'' are given by 
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FIG. 2. A combination of the 99% C.L. bound on |[/^4p [HI 
from CDHS and atmospheric neutrino data (left) and a 99% 
C.L. null result at Daya Bay can rule out the red points of the 
99% C.L. region favored by a joint analysis of LSND and Mini- 
BooNE antineutrino data in the 3+1 scenario i5l (right). The 
grey points survive the joint constraint, but not the ICARUS 
exclusion [9]. 



where S and R (=28,000/year) are the number of signal 
events from the source and the number of reactor back- 
ground events, respectively. 

The number of signal events (in all 4 detectors) with 
sterile neutrino oscillations is obtained by scaling the 
number of events for the case: 
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where An/ AE^s and An/ Ax are normalized event dis- 
tributions binned in visible energy and position, respec- 
tively, and Ntot = 38, 000 is the total number of events 
for the case in one year. The positron's energy in an 
inverse neutron /3-decay event is — (m„ — rrip). Subse- 
quent pair annihilation in the scintillator produces visible 
energy, 

Ey,s =E^- (to„ - mp) + me ~ - 0.8 MeV . (6) 

3+1. We checked that in the 3+1 scenario our pro- 
cedure yields a 95% confldence level (C.L.) sensitivity 
that is comparable to that of Ref. [7] for Smli < 2 eV^ . 
The oscillation amplitude that fits the global SBL data 
is given by 



Sin^ 20SBL =4|C/e4|'|C/,,4|' 



(7) 



(l + a,)5*7^" + (l + a,)i?. 



(4) 



Daya Bay data could push |C/e4| down far enough that 
the value of \U^4\ needed to obtain an amplitude that 
explains the SBL data could conflict with the current 
bound on | 1 shown in the left panel of Fig. [2] 

Since a meter-baseline measurement at Daya Bay will 
be independent of the earlier data, it is reasonable to 
impose the constraint on [7^4 as a prior. Then, Daya 
Bay can rule out most of the allowed region from a fit to 
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FIG. 3. 95% C.L. exclusion contours for the scenario, for 
\Ue4\ and lUebl above \U\min in steps of 0.01. The region 
above each contour is excluded. 
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\u,4 




\u,4 




KMS (3+2) [5] 


62 


0.47 


0.87 


0.128 


0.138 


0.165 


0.148 


KMS (1+3+1) [S] 


68 


-0.47 


0.87 


0.129 


0.142 


0.154 


0.163 


GL (3+2) g] 


78 


0.9 


1.61 


0.13 


0.13 


0.14 


0.078 


MM NH (3+2) [T0| 


64 


0.47 


0.87 


0.149 


0.127 


0.112 


0.127 


MM IH (3+2) 


80 


0.9 


1.61 


0.139 


0.122 


0.138 


0.107 


MMS (3+2) [n] 


55 


0.89 


1.76 


0.15 


0.07 


0.15 


0.15 



TABLE I. values for some global best-fit points (to data 
from SBL experiments), for a simulated dataset with no os- 
cillations at Daya Bay. 'MM NH/IH' is the 'minimal model 
(with normal/inverted mass hierarchy)' of Ref. [10] and 
'MMS' is the 'minimal seesaw model' of Ref. [11) . 



LSND and MiniBooNE antincutrino data; see the right 
panel of Fig. |2] 

Five neutrinos. We first consider Daya Bay's sen- 
sitivity to the scenario without recourse to spe- 
cific points, models or fits. We employ a grid in the 
{Sm1i,Sm'^i,\Ue4\,\Ue5\) parameter space, place a prior 
on the size of the mixing, min(|J7e4|, |J7e5|) = mm in 
steps of size 0.01 from 0.10 to 0.15, and suppose a null 
result at Daya Bay. The 95% C.L. sensitivity in the 
{617111,6171^1) plane is shown in Fig. [s] As mentioned 
before, P^e does not depend on the signs of the mass- 
squared differences. So the results of Fig. |3] apply to 
both the 3+2 and 1+3+1 spectra. 

We now specialize to models that are consistent 
with global neutrino data. In Table |lj we display Daya 
Bay's sensitivity to several best-fit points to SBL data 
in the case assuming that no oscillations are seen in 
the Daya Bay dataset. These points would be completely 
excluded by Daya Bay because of their sizable Ue4 and 

To examine Daya Bay's capability to probe the large biy 
parameter space, we use the globally allowed regions from 
an updated fit to the datasets listed in Ref. [5] in con- 
junction with data from the NOMAD [Hj and CDHS [TBI 
experiments |14j . The shaded areas of Fig. [4] are the glob- 




FIG. 4. The shaded regions are the 3a globally preferred re- 
gions in the {Srrili, Sm^i) plane in the 3+2 (left) and 1+3+1 
(right) scenarios ||14 . A null result at Daya Bay can rule 
out the regions above the blue solid and red dashed con- 
tours at the 95% C.L. and 3a, respectively. The exclusion 
regions are slightly different in the two panels because for 
each ((5m|i, (5m|i), the best-fit mixing matrix elements from 
the fit of Ref. [14] are difi'erent for the two mass spectra. 

ally allowed regions at Scr. We see that at least one 6m^ is 
close to 1 eV^ so as to explain the SBL data. All mixing 
parameters other than 6m/li and 6717^1 are marginalized 
over and assume their best-fit values. 

As the global fits favor significant — transitions, 
the mixing parameters tend to be large enough to be 
testable at Daya Bay. Figure |4] shows that Daya Bay can 
exclude the 3+2 and 1+3+1 scenarios as an explanation 
of the LSND/MiniBooNE anomaly at 3a. 

3+1 or 3+2? So far we have demonstrated that a 
null result at Daya Bay can significantly constrain sterile 
neutrinos. We now entertain the possibility that future 
data confirms their existence. Then, a pressing issue will 
be to ascertain whether the 3+1 or the 3+2 scenario is 
operative. Since scenarios with more eigenstates should 
be able to mimic those with fewer eigenstates, a good test 
of Daya Bay's discriminatory power is to fit 3+2 points 
to data simulated for 3+1. 

For every point on a grid in the {67n1i,6sBLi\Ue4:\) 
parameter space that lies within the 99% C.L. allowed 
region of the right panel of Fig. [2] and is also consistent 
with the 99% C.L. bound on \U^4\ in the left panel of 
Fig. [2J we simulate a dataset Nf^ . We then fit points in 
3+2 parameter space that are allowed at 3f7 (shown in 
the left panel of Fig. |4| to this dataset (using Eq. [3| , and 
find the 3+2 point with the minimum corresponding 
to that {67n1i,9sBL, \Ue4:\) point. For a given 6rnli, we 
repeat the procedure for other values of {9sbl, \ Ue4\) so 
as to find the global Xmin ^'^^ each (5m|j. Note that the 
best-fit 3+2 value of 6m1-^ need not be the same as the 
value for which 3+1 data was simulated. 

We plot Xmin versus SttiI^ in Fig. [5j The discrimi- 
nation between the 3+1 and 3+2 scenarios is better for 
small S7n1i. 
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FIG. 5. The degree to which Daya Bay can discriminate be- 
tween the 3+1 and 3+2 scenarios. We simulate an oscillation 
signal for points in the 99% C.L. region favored by LSND and 
MiniBooNE that are consistent with the 99% C.L. bound on 
|C/fi4| (see Fig. |2|, and fit the spectrum from points in the 
3(T region of the 3+2 parameter space (see the left panel of 
Fig. [4| to the simulated data. A 3(t discrimination is possible 
for <5mli < 0.5 eV^ 
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